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Abstract 

We study evolution of a flat Priedmann-Robertson Walker universe filled 
with a bulk viscous cosmological fluid in a higher derivative theory of gravity 
in the presence of time varying gravitational and cosmological constant. Cos- 
mological models admitting both power-law and exponential expansions are 
explored here in the presence of imperfect fluid described by full Israel and 
Stewart theory. We note some new and interesting cosmological solutions rel- 
evant for model building including present accelerating phase. In the case of 
power law, it is found that gravitational constant increases as the time evolves 
for a positive cosmological constant whereas it decreases for a negative cos- 
mological constant. The evolution of temperature of a viscous universe is also 
determined. 
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1 Introduction 

Recent astronomical observations of type la supernovae with redshift parameter Z < 
1 [1], Wilkinson Microwave Anisotropy Probe (WMAP) [2] etc. provided evidence 
that we may live in a low mass density universe {Q ~ 0.3) The predictions 
of the observations lead to a convincing belief in modern cosmology that a part 
of the universe is filled up with dark energy (f2 ~ 0.7), which may be addressed 
by a suitable cosmological constant. One of the recent predictions that the present 
universe accelerating is remarkable which is not properly understood yet. The present 
accelerating phase of the universe may be accommodated in cosmological models, 
obtained from theories : either by a modification of the Einstein's theory of gravity 
(GR) or by including a time varying cosmological constant. A number of literature 
appeared which addressed the present accelerating phase making use of exotic kind 
of fields in the matter sector of the Einstein's field equation. The study of late 
universe however, remains open to address the present issues in other theories. One 
of the early attempts to modify GR is based on adding curvature squared terms to 
the Einstein-Hilbert action which is known as generalized theory of gravity was used 
to obtain early inflation. The theory is found to have a number of good features 
important for understanding early universe. Starobinsky ^ shown that term 
in the Einstein Hilbert action admits inflation long before the advent of inflation 
is actually realized. However, the efficacy of the theory is known only after the 
seminal work on inflation by Guth [5j, who employed temperature dependent phase 
transition mechanism to obtain inflation. The inflationary scenario of the early 
universe is attractive and which may solve some of the outstanding problems in 
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cosmology [6]. It is known that higher order gravity with suitable counter terms viz., 
C^^P^C^ypS, B?, and cosmological constant (A) added to the Einstein-Hilbert action, 
one gets a perturbation theory which is well behaved, formally renormalizable and 
asymptotically free [TJ |8] . 

Cosmological models considering perfect fluid as a source of matter in the frame- 
work of higher derivative gravity have been studied in the literature [9] in order to 
obtain viable cosmological scenario of the early universe. In the absence of particle 
creation, particle number is conserved in the perfect fluid ( i.e., n!^ = 0, where n is 
particle number density ). It may be pointed out here that perfect fluid in equilib- 
rium generates neither entropy (S*." = 0) nor frictional heat flow as their dynamics 
is reversible. 

However, in the early universe a number of processes might have occurred leading 
to a deviation from perfect fluid assumption e.g. viscosity which is to be taken into 
account. It is known that real fluids behave irreversibly and therefore it is important 
to consider dissipative processes both in cosmology and in astrophysics. Some of the 
dissipative processes in the early universe responsible for viscosity may be due to 
the decoupling of neutrinos from the radiation era, the decoupling of matter from 
radiation during the recombination era, creation of superstrings in the quantum 
era, particle collisions involving gravitons, cosmological quantum particle creation 
processes and formation of galaxies [10]. It has been predicted from observations that 
a non negligible dissipative bulk stress on cosmological scales at the late universe 
phase might be important. The possible source of such viscosity may be due to 
(i) gaseous matter in the framework of relativistic gravity which may give rise to 
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internal self- interaction leading to a negative cosmic bulk pressure [TT] , (ii) deviation 
of the non relativistic particle in the substratum from dust. For a non-relativistic 
substratum cosmic anti-friction may generate a negative fluid bulk pressure which 
has been noted [12] in the framework of Einstein gravity. 

Since the pioneering work of Dirac [13], who proposed a theory with a time vary- 
ing gravitational coupling constant G{t), a number of literature [H] appeared with 
variable G{t) and A[t) in higher derivative theory of gravity in order to obtain cosmo- 
logical models accommodating the present cosmic acceleration. Recent applications 
of the apparent magnitude-redshift test, based on type la supernovae [21 [15], strongly 
favours a theory with a positive A term. A number of ansatze have been proposed for 
a dynamical decaying A with time [16]. Berman and collaborators [17\ pointed out 
that cosmological constant may vary as A ~ which seems to play a major role in 
cosmology. To obtain a time varying cosmological constant, one can assume different 
phenomenological relations proportional to the Hubble parameter, one such example 
[T7] considered in recent times is of the form A oc where H is Hubble parameter. 
Cosmological models with imperfect fluid have been explored in the presence of a 
dynamical cosmological constant and gravitational constant in Einstein gravity |18] . 
It is , therefore, important to look into viable cosmological models with variable G{t) 
and A{t) in higher derivative theory of gravity in the presence of imperfect fluid. 

To describe a relativistic theory of viscosity, Eckart [lOj made the first attempt. 
However, the theories of dissipation in Eckart formulation suffers from serious short- 
coming, viz., causahty and stability [20] regardless of the choice of equation of state. 
The problem arises due to first order nature of the theory, since it considers only 
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first order deviation from equilibrium. It lias been shown that the problems of the 
relativistic imperfect fluid may be resolved by including higher order deviation terms 
in the transport equation [21]. Israel and Stewart [22], and Pavon [23] developed a 
fully relativistic formulation of the theory taking into account second order devia- 
tion terms in the theory, which is termed as "transient" or "extended" irreversible 
thermodynamics (in short, EIT). The crucial difference between the standard Eckart 
and the extended Israel- Stewart transport equations is that the latter is a differential 
evolution equations, while the former is an algebraic relation. Extended irreversible 
thermodynamics takes its name from the fact that the set needed to describe non- 
equilibrium states is extended to include the dissipative variables (r, (). In irre- 
versible thermodynamics, the entropy is no longer conserved, but grows, according 
to the second law of thermodynamics (i.e., S^^ > 0). Bulk viscosity arises typically in 
mixtures either of different species or of the species but with different energies. The 
solutions of the full causal theory are well behaved for all the times. Therefore, the 
best currently available theory for analyzing dissipative processes in the universe is 
the full Israel-Stewart theory (FIS). Using the transport equation obtained from EIT, 
in addition to dynamical equation obtained from either Einstein gravity [TSl |25] or 
modified gravity [24J cosmological solutions are generally obtained. The motivation 
of this paper is to explore cosmological solutions in the modified theory of gravity in 
the presence of variable G{t) and A(t) with imperfect fluid described by FIS theory. 
It is interesting to study the behaviour of G{t) and A{t) in the presence of viscosity. 
The characteristic temperature of the universe may be determined in the presence of 
viscosity in the FIS theory. However, it may be pointed out here that in a number of 
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literature [2^ the evolution of temperature of the universe in the presence of viscosity 
from Gibbs equation is considered in the Einstein gravity. We also consider the above 
procedure in the framework of higher order gravity to determined the behaviour of 
temperature. 

The plan of this paper is as follows: in sec. 2, we give the gravitational action 
and set up the relevant field equations in the higher derivative theory of gravity. In 
sec. 3, cosmological solutions are presented. In sec. 4, distance modulus curves are 
presented. Finally, in sec. 5, we summarize the results obtained. 

2 Gravitational Action and Dynamical Equations: 

We consider a gravitational action with higher order term in the scalar curvature 
(R) containing a variable gravitational constant {G{t)) which is given by 

1 



I = - 



-fiR) + L„ 



-g d X [1] 



167rG(t)' 

where f{R) is a function of R and its higher power including a variable cosmological 
constant A(t), g is the determinant of the four dimensional metric and Lm represents 
the matter Lagrangian. 

Variation of the action (1) with respect to g^i, yields 

fniR) R,, - ^ f{R) g,, + fMR) (V^V.i? - g.uV^Vg^,) + 

fRRR{R) {Vf^RVuR - V'^RV^R g^u) = SnGit) T,,, (2) 

where is the covariant differential operator, fniR) represents the derivative of 
f{R) with respect to R and T^y is the effective energy momentum tensor for matter 



6 



determined by L^. We consider a flat Robertson- Walker spacetime given by the 
metric 

els'" = -df + a^t) [dr^ + r\d9^ + sin^9d(l)^)] (3) 

where a{t) is the scale factor of the universe. The scalar curvature for a flat universe 
is 

R^-6[H + 2H'^] (4) 

where if = ^ is the Hubble parameter and an overdot represents derivative with 
respect to cosmic time {t). The trace and (0,0) components of eq. (2) are given by 

RfniR) - 2f{R) + 3fRR{R) (^R + 3^^) + 3fRRR{R)R + 87rG{t) T = 0, (5) 

fRiR) Roo + biR) - ^fRR{R)-R + 87rG'(t) Too = 0. (6) 
^ a 

Let us consider a higher order gravity, namely, f{R) — R + aR^ — 2A(t). Using Eq. 
(3) in Eqs. (5) and (6), we get 



- 6a [2HH - + 6HH' 
and the conservation equation becomes 
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p+3(p + p)i/=-gp+^), (8) 

where p and p are the energy density and pressure of the perfect fluid respectively. 
Equations (7) and (8) are the key equations to study cosmological models with a 
perfect fluid in the presence of time varying G and A. To include, the effect of 
viscosity in the above, the perfect fluid pressure in eq. (8) is replaced by an effective 
pressure Peff, which is given by Peff = p + 11, where p is isotropic pressure and 11 



7 



is the bulk viscous stress. In EIT, the bulk viscous stress 11 satisfies a transport 
equation given by 



U + rtl^ -3CH - ^tU 



(9) 



where ( is the coefficient of bulk viscosity, r is the relaxation coefficient for transient 
bulk viscous effects and T > is the absolute temperature of the universe. The 
parameter e takes the value or 1. Here e = represents truncated Israel-Stewart 
theory and e = 1 represents full Israel-Stewart (FIS) causal theory. One recovers the 
non-causal Eckart theory for r — 0. The conservation eq. (8) including viscous ffuid 
is given by : 

^ + 3(p + p+n)i/ = -gp + ^). (10) 

For a constant G, A and 11 = 0, eq. (10) reduces to the usual continuity equation for 
a barotropic fluid. We consider an equation of state for the isotropic fluid pressure 
given by 

p=(7-l)p (11) 
where 7 (1 < 7 < 2) is a constant. The deceleration parameter (q) is related to H 



«-^f4)-l. (12) 



as 

dt \H. 

The deceleration parameter is negative for accelerating and positive for decelerating 
phase of the universe. The temperature of the universe is defined via the Gibbs 
equation which is given by 

rd5 = d(9+pd(i). (13) 
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The behaviour of temperature in the universe is obtained through Gibbs integrabihty 
condition 

dT . .dT dp 

7-1 

For a barotropic fluid the temperature follows a power-law which is T ~ p 7 . The 
above temperature may be determined using Gibbs integrabihty condition also 

(15) 



T 
T 



dp\ n dT' 



To determine the temperature of a viscous universe eq. (9) may be employed here, 
in addition to, Gibbs integrabihty condition (15). Both of them are considered here 
to investigate viable cosmological models. 

3 Cosmological Solutions : 

The system of eqs. (7), (9)-(ll) is employed to obtain cosmological solutions. The 
system of equations is not closed as it has eight unknowns (p, 7, r, G A, a{t), 
T) to be determined from four equations. We assume the following widely accepted 
adhoc relations 

c = /?p^ r = (16) 

where C>0,r>0,/3>0 and s > 0. We consider in this paper a known variation 
of A with Hubble parameter (H) oi the form, A = 3mH^ where m is dimensionless 
constant to be determined from the dynamical equations. In the next section we 
explore cosmologies with power law and exponential expansion respectively. 
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3.1 Power-law model: 

In this case we consider a power law expansion of the universe given by 

a{t) = aot^ (17) 

where oq and D are constants which are to be determined from the field equation. 
The accelerating mode of expansion (g < 0) of the universe is obtained for D > 1. In 
the absence of particle creation ( i.e., n!^ = ) eq. (10) may be decoupled as follows: 

p + 3-fpH + 3UH = 0, (18) 

Sir G p+ A ^0. (19) 
Using eqs. (7), (17) and (19) we obtain 



f + p2 , (20) 



where we replace p2 — ^^"'i^^^^ and Go — const. The initial value of G{t) is deter- 
mined in terms of coupling parameter (a), which however vanishes when D — ^. In 
this case the gravitational parameter G{t) increases with time for A > and it de- 
creases for A < 0. It is evident that the model admits a constant G when A vanishes. 
The energy density evolves as 



P = Po 



r + P2 r^ (21) 



where po — 3(1 — m)D^ (with SttGo = !)■ For physically realistic solution ( i.e., 
p > 0) the upper boundary on cosmological constant is A < (i.e., m < 1). 
For m — \ the variation of energy density is determined by term only for which 
it decreases as (p ~ t"''), whereas it is independent of a. The bulk viscous stress 
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obtained from eq. (18) is given by 



n = - 



Uot' + U2 t' + p2 '""'r^ (22) 



where Hq = D{3-fD{l - m) - 2) and U2 = lSaD{2D - l){3-fD - 4). For physically 

reahstic solution bulk viscous stress is essentially negative, which demands IIo > 

9 2 
i.e., D > ir-fr — T. We note that the bulk viscous stress decreases as 11 ~ t ^-^ 

' 37(1— m) I I 

when (i) q; = and D 7^ | or (ii) a 7^ and D — \. It is evident that for m < 1, 
bulk viscous stress (|n|) will decrease with time. One of the advantage of the FIS 
theory is that one can determine the evolution of the temperature in this case. For 
FIS theory we use e = 1, eq. (9) reduces to a differential equation given by 

^ = 3a-^^ + 5^ + p-i + ™ (23) 

T p n ^ ^ n ^ ' 

Using eqs. (17), (21) and (22) in eq. (23) we get temperature of the universe, which 
is 

T^T.^er-^^^e^h^-'^^ (24) 
P 

where To stands for a constant. The temperature of the universe have to go through 
Gibbs integrability condition (15) for a viable solution. We note the following: 

(i) For s — \^ m — \ and a = 0, we note power law decrease of the temperature, 

1 

which is T = To t""^ where a\ — + 4 — 3D — The decreasing mode of 
temperature is ensured for ol\ > 0. Putting the expression of temperature (T) in 
the presence of viscosity in Gibbs condition (15) one can obtain ai — Here 
the negative value of bulk viscous stress (11) is obtain for IIq > (i. e., 37D > 4). 
However, In the absence of viscosity the variation of temperature can be obtained 
from eq. (15) which yields T = Tot-^^^^-^\ It is evident from fig. (1) that the 
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temperature of the universe is more for a universe filled viscous fluid at a given time 
compared to that of a universe without viscosity. 

0.5 
0.4 
0.3 
0.2 
0.1 
0.0 

2 4 6 8 10 

t 

Figure 1: shows the plot of ^ vs t for a different value of 7 with D = |. Here the 
solid line represent the variation in the absence of viscosity while the dashing line 
represent that in the presence of viscosity. 

It is also evident that the temperature of the universe is higher for lower values 
of 7 at a given instant of time. 

(ii) For s = |,m=|,a7^0 and D ^ the temperature of the universe evolves 

1 1 

o ^ 3 

as T = Tot^^^"^" "2 +PP2>[Uot^ + U2r~^ ' [t^ + P2] ~^ . At the later 
stage of evolution of the universe (t >> 112 and t » P2) one can obtain the power 

law variation of temperature. Using Gibbs integrability condition the temperature 

-12(7-1) 

evolution of the universe becomes T = Trd 1 . Here the condition D > - ensure 

u ^ 

the negative value of bulk viscous stress. 

(iii) For s = the temperature of the universe in GR {a = 0) may be evaluated 
using Gibbs integrability condition (15) and eq. (9), we obtain the following variation 
of temperature : 

-2(7-1) , , 

T = Tot^^(^, (25) 
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which is decaying if7> l,m< lor7< l,m> 1. Here the negative value of 



bulk viscous stress (11) is obtain for Ho > ( i. e., D > 



37(1— m) 



). The condition 



D > also implies that in an expanding universe the temperature decrease 

less rapidly in the presence of viscosity compared to that when viscosity is absent. 
The plot of T Vs. t in fig. (2) shows that at a given time, a universe with higher 
temperature is possible when the cosmological constant is gradually decreased in 
GR. It permits a universe with late acceleration. 
Using the eq. (25), One can obtain the temperature of the viscous universe at the 

0.5 r 



0.3 - W 



0.1 - 



0.0 L 




m=0.3 



40 60 
t 



Figure 2: shows the plot of ^ Vs. t for a different value of cosmological constant 
(i.e., m) for 7 = |. 



different epoch of cosmological evolutions of the universe. To obtain the temperature 
at present epoch we set the initial condition Tq = 1 x 10^° at t = 1 [26] i.e., during the 
cosmological epoch of decoupling of neutrinos from the cosmic plasma. Considering 
age of the present universe ~ 10^° years ~ 3.15 x 10^^ second, 7 = | (radiation) 
and cosmological constant A = 3 x 0.085 x if^(z.e.,m = 0.085), eq. (25) leads to 
T ~ 2.74K which is in fair agreement with observed value T ~ 2.72K from CMBR. 
So, the solution admits observed value of temperature from CMBR in the presence 
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of a small value of cosmological constant with an late accelerating universe. 

3.2 Exponential models: 

The set of eqs. (7), (9), (18)-(19) also admit cosmological solution with a universe 
without singularity . We discuss two such cases. 

Case I : We note that Hubble parameter satisfying a differential equation 

H = nH-^-H\ (26) 



is permitted with rj = y-^- This admits a ever expanding universe with no sin- 
gularity, such a solution admits emergent universe scenario |27]. In this case the 
corresponding variation of gravitational constant and energy density are given by 

G = Goexp[bH-% (27) 

p = PoH^exp[-bH-% (28) 

where po = ^ (with SttGq = 1) and b = ^ and Gq = const. It is evident that the 
model admits a constant G when A = 0. The bulk viscous stress is obtained from 
eq. (18), which yields 

n = (2 - 7)p - -pH-^ \4:r]H^ + 2br] - SbH] . (29) 
3 

We note that without a cosmological constant (A = 0), one obtains a realistic solution 
when Hubble parameter satisfies an upper bound H < ^2~'r) ■ '^^^ temperature 
evolution of the universe in full causal theory is obtained from eq. (23), which yields 

T -n p-^(^) (-if]) 
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where f{H) = J ^^dt and g{H) = | / p^~^dt. The above expression of the tem- 
perature also have to satisfies the Gibbs integrabihty condition (15). We note the 
following : 

(i) For s = |, /3 = ^y^, in the absence of cosmological constant and for stiff fluid the 
temperature of the universe evolves as T = TqH'^{H — ^)~^. When H <^ the 
evolution of temperature becomes T ~ H^. We obtain here a decreasing mode of 
temperature for an exponentially expanding universe. 

The scale factor of the universe is obtained on integrating eq. (26) which is given 

by 

2 

a{t)± = ai ± 02 e''* ' . (31) 

The solution (a+) is important for building emergent universe scenario [27]. It has no 
singularity and the universe originated from a static state in the inflnite asymptotic 
past (t — )■ — oo). In this case, we note that the solution (a(t) = a(t)+) represents a 
universe which begins with a finite size in the past and grows exponentially. However, 
initially at t ~ 0, the universe is matter dominated which subsequently emerges to 
an accelerated phase of expansion. We note that the temperature of the universe 
decreases in this case. The decreasing mode of temperature obtained in the FIS 
theory may be relevant for the later evolution which will be taken up else where for 
a detail scenario. 

Case II : The dynamical equation admit de Sitter solution with H = The 
scale factor of the universe evolves as a(t) = aoe"*, for a sufficient infiation to solve 
cosmological problems, de Sitter phase should exit after an epoch At > We note 
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de Sitter phase with 

G = const., p = const., A = const. (32) 

and ( = const, and r = const, for m 7^ 1. In this case that bulk viscosity also remains 
constant throughout the inflationary phase. The temperature of the universe in this 
case is found to be a constant ( similar to that obtained by A. Beesham in Ref [18]) 



for (3 



7 p 



(2-7)r,- 

We also note that for H = ^r] and m = 1, one ends up with A = const., G = const., 
p = 0, C = 0, r = and T = 0, which is not physically relevant. 

4 Distance Modulus Curves : 

We now probe late universe with exponential and power law expansion taking into 
account observational results. The distance modulus is p = Slogd/, + 25, where the 
luminosity distance = ti(1 + z) a(to) and z represents the red shift parameter, 
where 1 + z = ^t^- We determine ri from 

a(ti) ^ 

(33) 



^ ti ^ ' 

2r}t 

At the late time, the scale factor becomes a{t) ~ a2e~. For exponential expansion 
of the flat universe, the distance modulus relation is given by 

M^)=51og(^y-^) + 25, (34) 

where Hq = y. For power law expansion (a(t) ~ t^) with k = 0, the distance 
modulus relation is given by 

f , D ^ 1 (1 + z) , , , D-i A , , 

M^) = 51og [ijfj''-^^ii^ + - 1) 1 + 25. (35) 
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The observed values of fi{z) at different z parameters [28] given in table 1 are em- 
ployed to draw the curves corresponding to the exponential and power law expansion 
of the universe discussed above. The plot are shown in figs. (3) and fig. (4) which 
matches with observations perfectly. 

48 
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0.0 0.5 1.0 1.5 

z 

Figure 3: shows the plot of /i Vs. z for supernova data (dashing line) and for the 
exponential expansion (solid line) with t] = 0.00038. 
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Figure 4: shows the plot of fi Vs. z for supernova data (dashing line) and for the 
power law expansion (solid line) with Hq = 10~^ and D = |. 
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z 


supernova /i 


Exponential /j, 


Power law fi 


0.038 


36.67 


35.961 


36.543 


0.014 


33.73 


33.742 


34.343 


0.026 


35.62 


35.112 


35.703 


0.036 


36.39 


35.839 


36.423 


0.040 


36.38 


36.077 


36.657 


0.050 


37.08 


36.582 


37.154 


0.063 


37.67 


37.110 


37.673 


0.079 


37.94 


37.634 


38.184 


0.088 


38.07 


37.887 


38.43 


0.101 


38.73 


38.212 


38.745 


0.160 


39.08 


39.324 


39.814 


0.240 


40.68 


40.349 


40.783 


0.300 


41.01 


40.936 


41.329 


0.380 


42.02 


41.579 


41.921 


0.430 


42.33 


41.925 


42.235 


0.490 


42.58 


42.298 


42.572 


0.526 


42.56 


42.504 


42.757 


0.581 


42.63 


42.797 


43.018 


0.657 


43.27 


43.165 


43.345 


0.740 


43.35 


43.530 


43.665 


0.778 


13.81 


13.686 


13.801 


0.828 


43.61 


43.881 


43.971 


0.886 


42.91 


44.096 


44.158 


0.949 


43.99 


44.316 


44.348 


0.970 


44.13 


44.388 


44.490 


1.056 


44.25 


44.665 


44.646 


1.190 


44.19 


45.061 


44.983 


1.755 


45.53 


46.403 


46.104 



Table 1: Red shift data from Supernova, power law model and exponential model. 
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5 Discussions : 

In this paper we study both power law and exponential behaviour of the universe with 
viscosity separately in a higher derivative theory of gravity considering a time varying 
cosmological and gravitational constant. We note that for a power law evolution of 
the universe, one obtains an increasing mode of gravitational constant with a positive 
A but a decreasing mode of gravitational constant results for a negative A. We note 
a physically reahstic (po > 0) solution for A < 3H'^. For A = the variation 
of energy density (p ~ t""^) is determined by term only. However, the energy 
density does not depend on the coupling parameter a in the gravitational action. 
We determine the characteristics temperature of a viscous universe from eq. (9), 
which also follows from Gibbs integrabihty condition. Fig. (1) shows the variation 
of temperature for different values of 7 in the presence or in the absence of viscosity 
in an accelerating universe (g > 0). It is evident that the higher value of 7 leads 
to a universe with lower temperature at a given instant of time. The evolution 
of temperature of a viscous universe is found to be more than that in a universe 
without viscosity. Fig. (2) shows the variation of temperature for different values 
of A (determined by m). It is evident from fig. (2) that the rate of decrease of 
temperature is higher for larger values of cosmological constant in the presence of 
viscosity. Here we obtain a interesting solution which suggests that the present 
temperature of the universe is T ~ 2.745K, which is in fair agreement with observed 
value T ~ 2.72K from CMBR, for a radiation dominated late accelerating universe 
in the presence of a positive cosmological (A = 3 x 0.085 x H'^) in GR. We also note 
cosmological solutions which admit a universe originating from singularity free state. 
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One of the solution corresponds to emergent universe [27] which is interesting. We 
obtain a de Sitter solution with A ^ 3H^. The figs. (3)-(4) show the plot of distance 
modulus (/i) vs red shift parameter {z). The figures indicate that power- law and 
exponential model support the present observational data perfectly well. 
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